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Abstract. This work presents a hybrid approach to solve the maximum 
stable set problem, using constraint and semidefinite programming. The 
approach consists of two steps: subproblem generation and subproblem 
solution. First we rank the variable domain values, based on the so- 
lution of a semidefinite relaxation. Using this ranking, we generate the 
most promising subproblems first, by exploring a search tree using a lim- 
ited discrepancy strategy. Then the subproblems are being solved using 
a constraint programming solver. To strengthen the semidefinite relax- 
ation, we propose to infer additional constraints from the discrepancy 
structure. Computational results show that the semidefinite relaxation 
is very informative, since solutions of good quality are found in the first 
subproblems, or optimality is proven immediately. 



1 Introduction 

This paper describes a hybrid method to solve the maximum weighted stable set 
problem, or stable set problem 1 in short. Given a graph with weighted vertices, 
the stable set problem is to find a subset of vertices of maximum weight, such 
that no two vertices in this subset are joined by an edge of the graph. In the 
unweighted case (when all weights are equal to 1), this problem amounts to the 
maximum cardinality stable set problem, which has been shown to be already 
NP-hard [21]. 

We propose a two-phase approach to solve the stable set problem, either 
with or without proving optimality. The first phase generates subproblems based 
upon a semidefinite relaxation, the second phase solves the subproblems using 
constraint programming. Concerning the first phase, given a model for the stable 
set problem, we solve its semidefinite relaxation. The solution provides us frac- 
tional values for the variables of the model. These fractional values are a good 
indication for the optimal (discrete) values of the variables. Hence we divide se- 
lected variable domains in two parts: a 'good' subdomain and a 'bad' subdomain. 
By branching on these subdomains using a limited discrepancy strategy [13], we 

1 Alternative names for the stable set problem are vertex packing, coclique or inde- 
pendent set problem. 



obtain first a very promising subproblem, and subsequently less promising sub- 
problems. 

The second phase consists of the solution of the subproblems. Since they 
are much smaller than the original problem, we can easily solve them using a 
constraint programming solver. 

As computational results will show, the semidefinite relaxation is quite in- 
formative. In several cases we can simply round the solution of the relaxation 
and obtain a provable optimal solution already. Otherwise, we are likely to find a 
good solution in one of the first subproblems. Using a limited number of subprob- 
lems to investigate, we yield an incomplete method producing good solutions. 
In order to obtain a complete search strategy, we need, in principle, to generate 
and solve all possible subproblems. A good upper bound is necessary to prove 
optimality earlier. For this reason we investigated the use of 'discrepancy cuts' 
that can be added to the semidefinite program to strengthen the relaxation and 
thus prune large parts of the search tree. However, computational results will 
show that they can not be applied efficiently on the instances we considered. 

The outline of the paper is as follows. The next section gives a motivation 
for the approach proposed in this work. Then, in Section 3 some preliminaries 
on semidefinite programming are given. In Section 4 we introduce the stable 
set problem, integer optimization formulations and a semidefinite relaxation. A 
description of our solution framework is given in Section 5. Section 6 presents 
the computational results. This is followed by an overview of related literature 
in Section 7. Finally, in Section 8 we conclude and discuss future directions. 

2 Motivation 

Let us first motivate why a semidefinite relaxation is used rather than a linear 
relaxation. Indeed, one could argue that linear programs are being solved much 
faster in general. However, for the stable set problem, linear relaxations are not 
very tight. Therefore one has to identify and add inequalities that strengthen 
the relaxation. But it is time consuming to identify such inequalities, and by 
enlarging the model the solution process may slow down. 

Several papers on approximation theory following [9] have shown the tight- 
ness of semidefinite relaxations. However, being tighter, semidefinite relaxations 
are more time-consuming to solve than linear relaxations in practice. Hence one 
has to trade strength for computation time. For large scale applications, semidef- 
inite relaxations can often be preferred as the relaxation of choice to be used in 
a branch and bound framework. Moreover, our intention is not to solve a relax- 
ation at every node of the search tree. Typically, we only solve a relaxation at 
the root node of the first phase (its solution is used to identify the subdomains), 
and optionally at the root node of each subproblem (in order to strengthen the 
upper bound). Therefore, we are willing to make the trade-off in favour of the 
semidefinite relaxation. 

Another point of view is the following. Although semidefinite programming 
has been developing for many years now in the operations research community, 



no efforts of integration or cooperation with constraint programming have been 
made to our knowledge. Application of semidefinite programming to problems 
typical to constraint programming, as was done in the papers on approximation 
algorithms mentioned in Section 7, is not yet hybrid problem solving. In this 
paper a first step is being made. The solution of the semidefinite relaxation is 
used to identify promising subdomains, and also produces a tight upper bound 
for the constraint programming solver. On the other hand, the solutions found by 
the constraint programming solver serve as a lower bound inside the semidefinite 
programming solver. 

3 Preliminaries on semidefinite programming 

Combinatorial optimization problems that are NP-hard are often solved with 
the use of a relaxation. Typically, the integral restriction on the decision vari- 
ables is relaxed to obtain a continuous relaxation, which can be polynomially 
solvable. Here we introduce semidefinite programming as an extension of the 
more common linear programming. Both paradigms can be used as a continuous 
relaxation. 

In linear programming, combinatorial optimization problems are modeled in 
the following way: 

max c T x 

s.t. Ax < b (1) 

x > 

Here x e K™ is an n-dimensional vector of decision variables and c£l"a cost 
vector of dimension n. The (m x n) matrix A G ffi> mxn and the m-dimensional 
vector b € R m define m linear constraints on x. In other words, this approach 
models problems using nonnegative vectors of variables. 

Semidefinite programming makes use of positive semidefinite matrices of vari- 
ables instead of nonnegative vectors. A matrix A e K" x ™ is said to be positive 
semidefinite (denoted by A > 0) when y T Xy > for all vectors y G K". Semidef- 
inite programs have the form 

max tr(CA) 

s.t. tr(A,A) = bj (j = 1, . . . , m) (2) 

x y o. 

Here tr(A) denotes the trace of A, which is the sum of its diagonal elements, 
i.e. tr(A) = Y^i=i^a- The cos ^ ma trix C € R nxn and the constraint matrices 
Aj e M" x " are supposed to be symmetric. The m reals bj and the m matrices 
Aj define again m constraints. 

We can view semidefinite programming as an extension of linear program- 
ming. Namely, when the matrices C and Aj (j = 1, ...,m) are all supposed 
to be diagonal matrices, the resulting semidefinite program is equal to a linear 
program. In particular, then a semidefinite programming constraint AjX = bj 
corresponds to the linear programming constraint defined by the j-th row of A 
and the j-th entry of b. 



Applied as a continuous relaxation, semidefinite programming in general pro- 
duces solutions that are much closer to the integral optimum than linear pro- 
gramming. Untuitively, this can be explained as follows. The problem is 'lifted' 
to a higher dimension, where one can exploit more of the structure of the prob- 
lem. Unfortunately, it is not a trivial task to catch this structure, and otain a 
good semidefinite program for a given problem. 

Theoretically, semidefinite programs have been proved to be polynomially 
solvable using the so-called ellipsoid method (see for instance [10]). In practice, 
nowadays fast 'interior point' methods are being used for this purpose (see [3] for 
an overview). Being a special case of semidefinite programming, linear programs 
are also polynomially solvable using an interior point method. However, they 
are often solved with a special linear programming solver, the Simplex method. 
Although this method can have an exponential behaviour in theory, in practice 
it is the fastest method known for solving linear programs. 

4 The stable set problem 

In this section, the stable set problem is formally defined, and formulated in 
two different ways as an integer optimization problem. From this, a semidefinite 
programming relaxation is inferred. 

4.1 Definition 

Consider an undirected weighted graph G = ( V, E), where V is the set of vertices 
and E the set of edges G V, i ^ j} of G. Let |V| = n and \E\ = m. To 

each vertex i G V a weight lu* € M is assigned. A stable set is a set S C V such 
that no two vertices in S are joined by an edge in E. The stable set problem is the 
problem of finding a stable set of maximum total weight in G. This value is called 
the stable set number of G and is denoted by a(G). The maximum cardinality 
(or unweighted) stable set problem can be obtained by taking all weights equal 
to 1. 

4.2 Integer optimization formulation 

Let us first consider an integer linear programming formulation. We introduce 
binary variables to indicate whether or not a vertex belongs to the stable set 
S. So, for n vertices, we have n integer variables a;, indexed by i € V, with 
initial domains {0, 1}. In this way, Xi = 1 if vertex i is in the stable set S, and 
Xi = otherwise. We can now state the objective function, being the sum of 
the weights of vertices that are in the stable set S, as Y^i=i WiXi - Finally, we 
define the constraints that restrict two adjacent vertices to be both inside S as 
Xi+Xj < 1, for all edges (i, j) <G E. Hence the integer linear programming model 
becomes: 

a(G) = max Yh=i w i x i 

s.t. Xi +Xj < 1 V(i, j) e E (3) 
x t g {0, 1} Vi G V. 



Another way of describing the same solution set is by the following integer 
quadratic program 

a(G) = max Ya=i w i x i 

s.t. XiXj = V(i,j) e E (4) 
xf = Xi \fi e V". 

Note that here the constraint iEj € {0, 1} is replaced by xf — x; t . This quadratic 
formulation will be used below to infer a scmidcfinitc programming relaxation 
of the stable set problem. 

In fact, both model (3) and (4) can be used as a constraint programming 
model. We have chosen the first model, since the quadratic constraints take more 
time to propagate than the linear constraints, while having the same pruning 
power. 



4.3 Semideflnite programming relaxation 

The integer quadratic program (4) gives rise to a semideflnite relaxation intro- 
duced by Lovasz [19], see Grotschel et al. [10] for a comprehensive treatment. 
The value of the objective function of this relaxation has been named the theta 
number of a graph G, indicated by Let us start again from model (4). As 

was indicated in Section 3, we want to transform the current model that uses 
a nonzero vector into a model that uses a positive semideflnite matrix to rep- 
resent our variables. In the current case, we can construct a matrix X G M™ x " 
by defining Xij — XiXj. Let us also construct a (n x n) cost matrix W with 
Wa = Wi for i e V and Wij = for all i ^ j. Since Xa = xf = x i: the objec- 
tive function becomes tr(WX). The edge constraints are easily transformed as 
XiXj = 0o X^ = 0. The first step in the transformation of model (4) can now 
be made: 

max tr(WX) 
s.t. X^ =0V(i,i) eE (5) 
xf = Xi Vi e V. 

This model is still a quadratic program, although reformulated. The problem 
remains how to model the last, very important, constraint. We need a mapping 
of the diagonal entries Xu — xf with the vector entries Xi. For this, we extend 
X with another row and column that contain vector x (both indexed by 0). For 
this reason, define the (n + 1) x (n + 1) matrix Y as 

/ 1 xi ■■■ x n \ 

Y - Xl 

X 

\x n J 

where the 1 in the leftmost corner of Y is needed to obtain positive semidefi- 
niteness. In this case we can express the required mapping as Ya = \Y i0 + \Y 0i 
(note that X and Y are symmetric), since then xf = Ya = jYo + ^Y M = a;,. 



The final step in the transformation consists of replacing the constraints on X 
by constraints on Y. In particular, instead of demanding X to be a product of 
nonnegative vectors, we restrict Y to be a positive semidefinite matrix. Namely, 
if the vector x represents a stable set, then the matrix Y is positive semidefi- 
nite. However, not all positive semidefinite Y matrices represent a stable set, in 
particular its values can take fractional values. 

In order to maintain equal dimension to 7, a row and a column indexed 
should be added to W, all entries of which containing value 0. Denote the 
resulting matrix by W. The theta number of a graph G can now be described as 

i?(G) = max tr(WY) 

s.t. Y u = iy i0 + 5*oi (i = 1, . . . , n) 
Y ij= V(i,j)€E 
Y±0. 

By construction, the diagonal value Yu serves as an indication for the value of 
variable x% (i € {1, . . . , n}) in a maximum stable set. In particular, this program 
is a relaxation for the stable set problem, i.e. i?(G) > a(G). Note that pro- 
gram (6) can easily be rewritten into the general form of program (2) . Namely, 
Yi — jYo + jYoi is equal to tr(AY) where the (n + 1) x (n + 1) matrix A 
consists of all zeroes, except for An = 1, A i0 = — \ and A^ = — 3, which makes 
the corresponding b entry equal to 0. Similarly for the edge constraints. 

The theta number also arises from other formulations, different from the 
above, see [10]. In our implementation we have used the formulation that has 
been shown to be computationally most efficient among those alternatives [11]. 
Let us introduce that particular formulation. Again, let x € {0, 1}" be a vector 
of binary variables representing a stable set. Define the n x n matrix X as 
X = -I— xx T . Furthermore, let the n x n cost matrix W be defined as 

f w i i= 3 
Wa = < . , 

for i, j e V. Then the following semidefinite program 

i?(G) = max tr(WX) 
s.t. tr{X) = 1 

Xn = ov(*,j) eE 

X>0 

gives exactly the theta number of a graph. Again, it is not difficult to rewrite 
this program into the general form of program (2). 

Program (7) uses matrices of dimension n and m + 1 constraints, while pro- 
gram (6) uses matrices of dimension n+1 and m + n constraints. This gives an 
indication why program (7) is computationally more efficient. 



(6) 



(7) 



Algorithm 1 Solution framework 
read problem 
set maximum discrepancy 

solve semidefinite program (7) — > upper bound 
round solution of (7) — > lower bound 
for i £ Vb do 

define Df ood and D^ d using solution of (7) 
end for 

set discrepancy = 

while lower bound < upper bound and discrepancy < maximum discrepancy do 
generate subproblem using LDS branching strategy on Df ood and D^ ad 
solve subproblem — > lower bound 
discrepancy = discrepancy + 1 

end while 



5 Solution framework 
5.1 Overview 

The two-phase solution approach proposed here is similar to the one described 
in [20]. In the first phase subproblems are generated, which are being solved in 
the second phase. A general overview of the method is presented in Algorithm 1. 

Let us first explain how we use the solution of the semidefinite program (7) to 
partition the domain Di of a variable Xi into Df°° and Z?^ ad (for i eV). As was 
stated before, the solution of program (7) assigns fractional values between and 
1 to the variables. Naturally, if for a variable Xi the corresponding fractional value 
Yu is close to 1, we regard 1 to be a good value for variable Xi. More specifically, 
we select the variable x% with the highest fractional value Yu, set Df ood = {1} 
and D\ ad = {0}, and mark it as handled. Then we mark all its neighbours j 
(with e E) as being handled, keeping their original domain Dj = {0,1}. 
This procedure is repeated until all variables are handled. For convenience, we 
partition V into two distinct sets V = {i e V \Df ood = {1}} and Vi = V \ V . 
Here V\ represents the set of neighbours of Vo, with Dj — {0, 1}. 

In a similar way, we can use the solution of the semidefinite relaxation to 
compute a first feasible integer solution. Namely, follow the same procedure, but 
now instantiate the selected variable Xj = 1 and set its neighbours Xj — 0. The 
objective value of this feasible integer solution is in many cases already equal 
to the (downward rounded) solution of the semidefinite relaxation. In that case, 
we have found an optimal solution and finish. In other cases, we can still use 
this first solution as a lower bound to be applied during the solution of the 
subproblem. 

Next, we explain how to generate subproblems using these subdomains. The 
generation of subproblems makes use of a tree structure of depth |Vo| in which 
we branch on Df ood versus D^ ad . The tree is traversed using a limited discrep- 
ancy strategy (LDS) [13]. LDS visits the nodes of a search tree differently from 
depth-first search. It tries to follow a given suggestion as good as it can. Branches 



opposite to the suggestion are regarded as discrepancies and are gradually al- 
lowed to be traversed. The first 'run' of LDS doesn't allow any discrepancies, 
the second allows one, and so on. This means that for a particular discrepancy 
k, a path from the root to a leaf is allowed to consist of maximally k right 
branches. Typically this method is applied until a limited number of discrepan- 
cies is reached (say 2 or 3), which yields an incomplete search strategy. In order 
to be complete, one has to visit all nodes, up to discrepancy d for a binary tree 
of depth d. 

In our case, the suggestion that should be followed are branches of the kind 
Df ood , while branches Z)j? ad are regarded as discrepancies. Hence, our first sub- 
problem is the subproblem defined by program (3), and with Xl £ Df od for all 
i £ Vq. The next \Vq\ subproblems have all Xi £ Df°° , except for one Xk £ D h £ d 
(i, k £ Vo). The next discrepancy generates ^(|Vb| 2 + \Vq\) subproblems, each of 
which contains two variables £ D^ d and Xk 2 £ D b j? d (ki,k 2 £ Vo), and so on. 
Since we expect to obtain a very good solution already in the first subproblems, 
we will only generate subproblems up to a certain maximum discrepancy. In our 
experiments the maximum discrepancy is chosen 2 and 4 respectively. 

Note that the first subproblem, corresponding to discrepancy 0, only con- 
tains one solution, namely the one that we obtain in our rounding procedure. 
By propagation of the edge constraints, all variables Xj £ Dj = {0, 1} are instan- 
tiated automatically to 0. Hence, the subproblem corresponding to discrepancy 
is obsolete in our current implementation. 



5.2 Adding discrepancy cuts 



In the case one needs to generate and solve subproblems up to a large discrep- 
ancy, it is preferable to proof possible suboptimality of a subproblem before 
entering it. Especially when the subproblems are still relatively large. This can 
be done in several ways. 

First, before entering a subproblem, we can identify variables which have a 
subdomain of size 1, namely those with i £ Vo. For those variables, one can 
add an extra constraint to the semidefinite program (7), being either X{ = 1 or 
Xi = 0. Then the semidefinite program can be solved again, and will in general 
provide a tighter bound, hopefully lower than the current lowerbound, in which 
case we have proven suboptimality. However, solving the semidefinite program 
each time before entering a subproblem is very time-consuming and this method 
will not be very practical. 

A better alternative would be to add a specific constraint, a discrepancy cut, 
that is valid for all subproblems of a given discrepancy. Recall that Vq = {i £ 
V |_Df ood = {1}}. Hence, all subproblems of discrepancy k consist of k variables 
Xi with Xi = 0, and \Vq\ — k variables Xi with Xi = 1 (i £ Vq). This gives rise to 



two discrepancy cuts, given discrepancy k: 
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We implemented both of them, and cut (9) gives the best results. Stated in terms 
of semidefinite program (7), the discrepancy cut looks like tr(AX) — \Vo\ — k, 
with Aij = 1 if cither i € Vo or j G Vq or both, and Ay = otherwise. As 
mentioned before, solving a semidefinite relaxation is relatively expensive, and 
one should make a tradeoff between its computation time and the gain in time 
of not solving the subproblcm. For the instances we considered, the time needed 
to solve a semidefinite program is always larger than the time needed to solve 
all subdomains we would like to proof suboptimal. However, these cuts might be 
helpful for larger instances. 

6 Computational results 

Our experiments are being done on a Pentium 1GHz processor, with 256 Mb 
RAM, using Windows XP. As constraint programming solver we use the ILOG 
solver library, version 5.1 [16]. As semidefinite programming solver, we use CSDP 
version 4.1 [5]. 

The first instances we consider are randomly generated weighted graphs with 
n vertices and m edges. The vertex weights are randomly chosen from a range 
of 1 up to n. The edge density is chosen such that the constraint programming 
solver has difficulties solving them. Namely, the more edge constraints we have, 
the easier it is solved by constraint programming. The name of the instances 
represent the number of vertices and the edge density, i.e. g75d015 is a graph 
on 75 vertices with an approximate edge density of 0.15. For these graphs, we 
have chosen to generate subproblems up to a maximum discrepancy of 4, based 
upon earlier experience. 

We also considered structured instances, obtained from problems arising in 
coding theory [23] . These are unweighted graphs, therefore we have set all weights 
equal to 1. For these graphs, we generate subproblems up to a discrepancy of 2. 

The results of our experiments are given in Table 1. It consists of three 
parts: the first part describes the instances, the next part gives the results of 
our approach (sdp and cp), the last part concerns the results of a sole constraint 
programming approach (cp alone). 

The columns in this table represent the following. An instance name has n 
vertices and m edges. For the part on our approach, the value of the semidefinite 
relaxation is the rounded solution of the semidefinite relaxation has value 
round, and best is the value of the best solution found. This best solution is 
found in a subproblem generated during discrepancy best discr. Note that we 
generate subproblems up to discrepancy 4 in all cases, as was mentioned in 
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Section 5. The time spent on solving the semidefinite relaxation is denoted by 
time sdp. The time spent on solving all generated subproblems is denoted by 
time subp. These values together form the total time. All times are measured in 
seconds. The number of all backtracking steps made during the search in our 
approach is collected in backtracks. Concerning the sole constraint programming 
approach, we report the best solution found (best), the total time spent during 
search, and the total number of backtracks. Note that we have set a time limits 
for the constraint programming solver, to create a fair comparison with our 
approach. They are 100 seconds for g50d005 up to gl50d010, 190 seconds for 
gl50d015 and 324 seconds for the structured instances. Best found solutions 
that are proven to be optimal are indicated by an asterisk (*). 

For the instances in Table 1 we only solved one semidefinite relaxation per 
problem, namely at the root node. The reason for this is that the time spent 
during the subproblem search is less than the time spent on computing another 
relaxation, as reported in the table. Therefore, we cannot gain time by adding 
discrepancy cuts and computing another semidefinite relaxation. 

In general, our method produces better bounds than the constraint program- 
ming approach alone. In many cases, the rounded solution of the semidefininitc 
relaxation is already optimal. However, note that there are two instances for 
which the constraint programming approach gives better bounds. Note also that 
the structured instances are handled quite well by our approach, while the con- 
straint programming approach produces very low quality solutions for the larger 
instances. 

7 Related literature 

Since the stable set problem is NP-hard, no complete (or exact) algorithm is 
known that solves the stable set problem in polynomial time. Many other tech- 
niques have been proposed, including approximation algorithms, heuristics, or 
branch and bound structured methods. A survey of different formulations, com- 
plete methods and heuristics for the maximum clique problem 2 is given by Parda- 
los and Xue [22] and, more recently, by Bomze et al. [4]. 

Although semidefinite programs can be solved in polynomial time theoret- 
ically, it lasted until a few years ago until fast solvers for this purpose were 
implemented. Until then, application inside a branch and bound framework was 
unrealistic. Still, solving a semidefinite program takes relatively much time, com- 
pared to solving a linear program. However, since semidefinite programming 
solvers are getting faster, semidefinite relaxations become a serious candidate to 
be used within a branch and bound framework, see for instance the paper by 
Karisch et al. [17]. 

A large number of references to papers concerning semidefinite programming 
are on the web pages of Helmberg [14] and Alizadeh [2]. A general introduction 

2 The maximum weighted stable set problem of a graph is equivalent to the maximum 
weighted clique problem of its complement graph. 



on semidefinitc programming applied to combinatorial optimization is given by 
Goemans and Rendl [8]. 

Another area that made semidefinitc programming useful in practice is that of 
approximation algorithms. In this field one tries to give a performance guarantee 
for an algorithm on a particular problem. In particular, the paper [9] by Goemans 
and Williamson uses a semidefinite relaxation and randomized rounding to prove 
such a performance guarantee for the maximum cut problem of a graph and 
satisfiability problems. Following this result numerous papers appeared, also 
concerning the approximation of satisfiability problems, including [12] and [18]. 

The solution structure of the current work, namely problem decomposition by 
branching on promising subdomains, is similar to the method described in [20]. 
In that paper a linear relaxation is used to identify promising values. Moreover, 
by exploiting the discrepancy structure of the method combined with reduced 
costs, suboptimality of subproblems can be proved very fast. 

Another hybrid approach using linear programming and constraint program- 
ming has been investigated by Ajili et al. [1] and El Sakkout et al. [6]. A subset 
of constraints is relaxed as a linear program in such a way that its solution is al- 
ways integral. The solution to the relaxation serves as a suggestion (a 'probe') for 
solving the complete program using a constraint programming solver. A probe 
is used to detect infeasibility, to remove inconsistent domain values and to guide 
the search. During search, many probing steps are being made. This results in a 
tight cooperation of the linear programming and constraint programming solver. 



8 Conclusions and further research 

We introduced a method that combines semidefinite programming and constraint 
programming to solve the stable set problem. Our experiments show that con- 
straint programming can indeed benefit greatly from semidefinitc programming. 
On instances that were very difficult to handle for a constraint programming 
solver, our hybrid method obtained very good results. 

The discrepancy cuts we proposed to strengthen the semidefinite relaxation 
could not be applied efficiently to the instances we considered. However, for 
larger instances they could be helpful. 

Further research in this direction would for instance be to obtain a filtering 
mechanism similar to the cost-based domain filtering for linear relaxations [7]. 
In [15], Helmberg describes such a procedure, called variable fixing, for semidef- 
inite relaxations. It would be interesting to see how his method can be applied 
in a constraint programming framework. 

Also, one could consider a different way of selecting promising values from the 
solution of the semidefinite relaxation. A strategy that incorporates randomized 
rounding possibly yields better results. This thought is motivated by the use of 
randomized rounding of semidefinite relaxations in approximation algorithms, 
as discussed in Section 7. 



Finally, this work has much in common with our previous work [20] . The un- 
derlying general principle of decomposing a problem into promising subproblems 
according to a certain heuristic is currently under research. 
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